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Abstract 

q^"' BES II data for the channel J /\P — ► 7(7r + 7r~7r + 7r~) are analysed from a sample of 58M 

J/4' hadronic interactions. They reveal further detail in the J PC = h sector for the 
4-7T channel. A broad — > component extending from 1500 to 2400 MeV is fitted. 
A narrower structure reveals a new J p = 0~ contribution requiring a resonance of mass 
1970 ± 25(stat) ± 60(syst) MeV and width T = 210 ± 25(stat) ± 60(syst) MeV. There is 
also tentative evidence for a 0~ resonance with mass 1560 ± (25) ±60(syst) MeV and width 
~ 280 MeV. In the + sector, there are clear contributions from /o(1500) and /o(1790). The 
observed /o(1790) is visibly distinct from /o(1710) and confirms the existence of /o(1790); 
its mass 1800 ± 25(stat) ± 20(syst) MeV and width 230 ± 30(stat) ± 25(syst) MeV are 
QJ ' better determined here than elsewhere. 
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The analysis reported here was done in the years 2001-3 in collaboration with three senior 
members of the BES collaboration. The work was funded by the Royal Society and Queen 
Mary College via a contractual agreement with the Chinese Academy of Sciences and the BES 
collaboration. This contract guaranteed access to BES 2 data for the purpose of publications. 
Since then, the BES management has refused to appoint internal referees, hence blocking a 
collaboration publication. The data illuminate the spectroscopy of J p = 0~ mesons and should 
^ ■ be in the public domain. There is a responsibility to publish work which has been supported 
from public funds with approximately £60, 000 . The wording of the present article is close to 
a document submitted late in 2008 to the BES management, but some further detail is added 
in Section 8. 

There have been several earlier studies of J/ty — > 7(7r + 7r~7r + 7r~). The Mark III [1] and DM2 
collaborations [2] observed a large 0~ signal and peaks at 1500, 1760 and 2100 MeV. Later, Expt. 
E760 at Fermilab found three peaks in rji] [3] with masses and widths very similar to those in 
J/^> —>■ 747r. In the rjr] channel, quantum numbers J p = 0~ are forbidden by Bose statistics. A 
re-analysis of the Mark III data [4] led to the conclusion that the peaks have quantum numbers 
J PC = ++ , but sit on a large h signal extending over the whole mass range up to 2400 MeV. 
BES I data led to similar conclusions [5], but added the identification of a broad 2 ++ signal at 
~ 2 GeV, decaying to /2(1270)cr and pp. 
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Here statistics a factor 7 higher are presented from BES II, using 58M J /iff hadronic inter- 
actions. The new data reveal further detail in the 0~ sector. There is clear evidence for a CT 
resonance at 1970 MeV. This new feature alters details of what may be fitted in this mass range 
to / (2100) and the J p = 2 + amplitude. There is also a strong pp peak at ~ 1600 MeV. Two 
alternative ways of fitting this peak are presented. The first alternative is that it may simply 
be a threshold pp cusp. The second investigates evidence that the peak is actually resonant. 
The distinction between the two alternatives depends rather delicately on the phase of the 0~ 
amplitude. In the J p = + sector, there is clear evidence for the existence of an / (1790) which 
is definitely distinct from / (1710). 

The layout of the paper is as follows. Section 2 concerns technicalities of data selection and 
reduction of experimental backgrounds. There is a troublesome background at high 4-ir mass 
from 7r°7r + 7r~7r + 7r~, mostly O2(1320)p and ai(1260)p. This background has little consequence 
below 2 GeV, but does obstruct separation of J p = 2 + and + from 2 to 2.3 GeV. 

Sections 3 and 4 present a fit to the An mass spectrum in slices of mass. The key point is that 
J /ty radiative decays are pointlike, because of the large cc mass and because photons interact at 
a point. The consequence is that only the lowest angular momentum states contribute, making 
the selection of partial waves clean, see Fig. 3. The dominant J p = 0~ channel, produced by an 
Ml transition, varies in intensity as P^, where P 7 is photon momentum. Superimposed on this 
broad signal is a narrow structure requiring the phase variation of a resonance at ~ 1970 MeV. 

For J p = + , there is a strong clear / (1790) peak, definitely distinct from f (1710), see 
Figs. 3(b) and (c). The 2 + channel is fairly weak. There is a definite signal in the mass range 
2-2.3 GeV, but it is difficult to say whether it is a single broad component or may be due to 
more than one resonance. The experimental background is the limiting feature in this mass 
range. Section 5 shows detailed fits to Air mass projections. 

Section 6 concerns the question of how to fit the dominant very broad 0~ component. From 
other data, it is known that there are threshold peaks, see Fig. 10, in uu, KK*, K*K* and 
00, as well as in the strong pp signal observed here. One or more of these peaks could be 
resonant, but it seems unlikely there is one resonance per channel. The alternative is that they 
are threshold cusps. The formula fitting the broad 0~ component is based on this conservative 
approach. Section 7 concerns fitted branching fractions. 

The fit to the 4tt mass spectrum from 1500 to 1650 MeV is not perfect in this approach. This 
mass range can be fitted significantly better if the strong peak observed in pp (and elsewhere in 
r/nn) is resonant. However, in the absence of strong interferences with other partial waves, this 
evidence for a resonance is tentative, and is discussed in Section 8. 

There is one feature of the present data which is in strong disagreement to an earlier pub- 
lication claiming an ^(1760) resonance which completely dominates data on J/$ ^ ^ujuj. It is 
suggested in Section 9 that an alternative way of fitting those data is in terms of a simple cusp 
due to the uu P-wave threshold. Section 10 summarises conclusions and the Appendix outlines 
the SU2 relation concerning Section 9. 

2 Event Selection and Backgrounds 

The data presented here were taken with the BES II detector. Full details of the detector and its 
upgrade are reported by Bai et al. [6], [7]. It has cylindrical symmetry around the intersecting 
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e + e~ beams. Its essential features are (i) a Main Drift Chamber (MDC) for the measurement of 
charged particles, (ii) time of flight detectors with a a of 180 ps, and (iii) a 12 radiation length 
Barrel Shower Counter comprised of gas proportional tubes interleaved with lead sheets. The 
MDC measures dE/dx. Together with the time of flight detectors, it separates tt and K up 
to ~ 700 MeV/c. Outside the 7 detectors is a magnet providing a field of 0.4T. This magnet 
is instrumented with muon detectors, but for present work they serve simply to reject fi + fi~ 
pairs. The Main Drift Chamber provides full coverage of charged particles for lab angles 9 with 
|cos#| < 0.84. 

Events are selected initially by demanding four charged tracks in the Main Drift Chamber 
with balancing charges, plus at least one photon in the lead-scintillator detectors. Any charged 
track with either time-of-flight or dE/dx favouring kaon identification is rejected. Charged 
tracks are required to have a transverse momentum > 60 MeV/c; this reduces background from 
J/^f — > 7r°7r + 7r~7r + 7r~. A geometrical cut is applied to ensure that all particles lie inside a region 
well inside the known detector geometry, so as to eliminate edge effects in detector performance 
and acceptance. The vertex is required to lie within a cylindrical region of radius 2 cm and ±20 
cm in length, with respect to the known intersection point of the e + e~ beams. 

Events have then been fitted kinematically to 7(7r + 7r~7r + 7r~), using all photon showers in 
turn; some photons may originate from interactions of charged particles in the detector, but any 
photons within a cone of 6° around a charged particle are rejected. Events with a x 2 probability 
> 5% are selected. It is also required that x 2 is better than that for likely background channels: 
(a) 47r ignoring photons, (b) 2747T if there is more than one photon and M(27) is within 50 
MeV of the 7r° mass. A small detail is that background from J/\I/ — > ujtt + tt~ has been rejected 
by a cut \M n +^- % o — > 25 MeV, where the n° is associated with missing momentum. Also, 
to remove a small background from J/\l> — > yK s K s , a double cut is made on the two 7r + 7r~ 
masses: \M w+7r -\ > 25 MeV. 

The Monte Carlo generates 300K events to simulate the BES detector and its acceptance; 
it uses the standard SIMBES package released in February 2002, in its version as of March 
2003. During discussions with the BES management, they have claimed that alterations to the 
calibration of the time-of-flight system after March 2003 could affect the selection of events. 
I have investigated this claim thoroughly and it has no foundation. Firstly, there is excellent 
agreement between dE/dx and time-of-flight measurements. Secondly, a search has been made 
for contamination of data due to misidentified / yK + K~tt + tt~ events. These are rather well 
separated from •yin events by kinematics as well as dE/dx and time-of-flight. For 60% of 
events, the kinematic fit alone fails to fit 'yK + K~-K + n~ by a difference of x 2 > 50. This is 
because production angles are sensitive to the different masses of kaons and pions. 

Fig. 1(a) shows the mass distribution for n^n^ pairs fitted to 747r with momenta P± > P3 
and P2 > P4. (Results for other 7r + 7r~ combinations are similar). There is an obvious p(770) 
signal; a small dip at 500 MeV arises from the cut against Kg. Next, these events have been 
fitted to all four combinations of ^K + K~it + it~ after dE/dx and time-of-flight cuts. Events 
fitting with a x 2 < 50 have the Kn mass distribution shown in Fig. 1(b). It is known from 
BES 1 data [8] that the / -/K + K~tt + tt~ channel contains a very strong i^*(890) signal. There 
is no evidence for this in Fig. 1(b), demonstrating that any contamination from / -fK + K~7T + 7T~ 
is < 1% of good 747T events. Likewise there is no structure in the KK mass spectrum of Fig. 
1(d). Fig. 1(c) shows the irir mass spectrum from events surviving the dE/dx and time-of-flight 
cuts but fitted as 7 K + K~ "7r + ~ir~ . One sees that the effect of the kinematic fit to the wrong 



3 




Figure 1: Mass projections for (a) tttt pairs from events fitted to 77r + 7r~7r + 7r~ , taking 
P(^i) > P(nt) an( i P( n 2 ) > P( n t)'i other panels refer to events fitted to all combinations 
of r yK + K~'K + 'K~ : (b) both Kn pairs, (c) 7r + 7r~ pairs, (d) KK pairs. 



K~it + it~ hypothesis is to shift the p(770) peak down in mass. 

Other background processes and the acceptance for data are studied using the standard 
SIMBES Monte Carlo simulation of the detector. The main background arises from the 7r°7r + 7r~7r + 7r 
channel after the loss of one photon. This is well known from all earlier publications. The 
7r°7r + 7r~7r + 7r~ channel has been studied by DM2 [9]. It contains a strong a 2 (1320)p signal, a less 
prominent ai(1260)p signal and thirdly a diffuse signal over all of phase space. The Monte Carlo 
has been used to generate a 2 (1320)p and ai(1260)p events with the acceptance of the BES 2 
detector. Background from these channels in ^Att data peak strongly at high An mass, because 
of the loss of low energy photons associated with low energy 7T°. This background dominates 
the 4-7T mass range > 2400 MeV/c 2 , allowing a firm determination of its magnitude. 

Both backgrounds produce clear p(770) — > tt°tt ± , see Fig. 2(a). In these events, the back- 
ground due to a 2 (1320)p leads to an a 2 (1320) peak in the pn mass distribution, shown in 
Fig. 2(b). Its magnitude is fitted freely and is consistent within errors with the branching 
fraction quoted by the Particle Data Group (PDG) [10]. The a± is sufficiently broad that the 
ai(1260)p background after cuts is consistent within errors with 5n phase space and is fitted 
like that. This background contribution is somewhat larger than DM2 quote, suggesting that 
further diffuse background close to 5ir phase space is also present. The combination of ai(1260)p 
and diffuse background is again fitted freely to the 4-7T mass range above 2.4 GeV. The sum of 
experimental backgrounds from a 2 (1320)p and a!(1260)p is shown below by the dotted curve in 
Fig. 5(d). Fig. 2(c) shows the observed tt°tt + tt~ mass spectrum in events fitted to it it + it~it + -7t~ . 
The uj peak disappears after cuts used to select ^An. 

Fig. 2(d) shows the whole mass distribution for it + it~it + it~ after the full event selection. 



4 




M(n°n + ri) (GeV) M(n + ri% + K) (GeV) 

Figure 2: Mass projections: (a)-(c) refer to events where charged particles alone are consistent 
with the 7r°7r + 7r~7r + 7r~ _ hypothesis, (a) M(7r°7r ± ), (b) M(3vr) after selecting M(7r°7r ± ) from 670 to 
870 MeV, (c) M(7r°7r + 7r~), all events; (d) M( y n + n^n + n^) after the full selection of 77T + 7r~7r + 7r~ 
events. 



Peaks are visible below a mass of 2.4 GeV. Above that, background dominates, so the analysis 
reported here is confined to the mass range below 2400 MeV. The total background up to 2.4 
GeV is 27.3% of the whole ^Ati sample. 

From Monte Carlo studies, estimated backgrounds from K + K~7T + 7T~ and n° K + K~ n + -k~ 
are < 2%; the latter background is peaked at high K + K~7t + tt~ masses, like the 7r 7r + 7r~7r + 7r~ 
background, and is mostly removed by selecting events with Att mass < 2400 MeV. 

3 Introduction to formulae for partial waves 

Required channels are pp, oa and /2(1270)cr. There are several parametrisations of the a 
amplitude with increasing complexity as its couplings to KK, r\r\ and 4-7T have been unravelled 
[11], [12], [13]. Only the coupling to tttt is important here. All have been tried in the analysis, 
but there is little to choose between them. The first formula parametrises tttt data adequately 
and gives slightly the best fit; it is the simplest and is used here. 

There might be contributions from resonances decaying to a 2 (1320)7r and/or ai(1260)7r. 
There is no evidence for from this source in mass projections of p~n after background sub- 
traction. Decays of 0~ to ai7r are forbidden by parity conservation. Adding decays of fo and / 2 
to ai7r from each resonance one by one gives negligible improvements in log likelihood (< 3). 

The angular momentum in decays of resonances will be denoted by L and the angular 
momentum in the production process J/ty — > 'jX will be denoted by V . The dominant process 
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Figure 3: J p = , + and 2 + components from a fit in 50 MeV slices of 4-n mass from 1350 to 
2350 MeV. (b) shows the sum of era and pp signals; the former dominates strongly. 



for production of X with J p = 0~ requires V — 1. The coupling of the photon to cc is a 
point interaction. Because the charmed quarks are massive, the interaction of gluons with cc 
is expected to be of very short range, with a radius parameter of order l/m c ~ 0.1 fm. The 
intensity for J/ty — > 70~ is then proportional to P 3 , (where P is the momentum of the photon 
and hence the 4n system), multiplied by a very slowly varying form factor discussed below. 
The data are indeed close to the predicted P 3 dependence. This has an important consequence: 
production with V > 1 should be strongly inhibited. Higher V values were tried in the fit, and 
in all cases are negligible. This leads to a major simplification in the analysis. 

For production of J p = 0~, the V — 1 operator is represented by the total momentum Kx 
of X in the J/^f rest frame. When it decays to PaPb, followed by pa — > ^1^2 an d Pb — > ^t 71 ! 
the vectors describing pa and pb are f^i = k\ — k 2 and K B = k 3 — /c 4 , where are momenta of 
pions in the rest frame of each p. [The 14 and 23 combinations also need to be added coherently 
in the program]. The full formalism for Lorentz invariant tensors is presented in Ref. [14]. It 
involves a Lorentz transformation of K\ and Kb to the centre of mass of X. Apart from this 
complication, the formula for production and decay of 0~ is given by K A A K B .K X - This form 
is highly distinctive, particularly the vital dependence on K X - 

The states with J p = + and 2 + are produced with V — 0. The + amplitude for 00 decays 
is isotropic; for decays to pp it is proportional to K^.Kb- The 2 + amplitude is given by the 
well-known spin 2 tensor: K\K^ — (1/3)(Ka-Kb)- All three J p components are cleanly and 
simply identified. There is no evidence for components with other J p , e.g. 1 + and 2~. A check 
on all amplitudes in the program is that they are orthogonal when integrated over 47r solid angle. 
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Figure 4: Differences in the slice fit between phases of (a) 2 + (pp) and + , (b)2 + [/ 2 (1270)o-] and 
2 + (pp). 



4 The Slice Analysis 

The data are fitted in 50 MeV bins of An mass and results are shown in Fig. 3. The + amplitude 
peaks near 1790 MeV and has a shoulder near 1500 MeV. For J p = 2 + , decays to both pp and 
/2(1270)er are observed. There is some cross-talk between these two final states, but both are 
consistent with a broad peak at ~ 2150 MeV, plus a small / 2 (1270) — > pp at masses 1270-1400 
MeV. The latter is weak and not well defined by present data. Its magnitude is known from the 
analysis of J/^f — > 77r + 7r~ and the 3.3% branching fraction of / 2 (1270) to p°p°, and is fixed to 
this value in the full fit described below; its phase is fitted freely. 

From Fig. 3(a), the J p = 0~ amplitude appears to be made of two components: peaks at 
1600 and 2100 MeV. However, the data cannot be fitted successfully like that. The reason is that 
there is clear evidence for phase variation consistent with a resonance near 1970 MeV interfering 
with the broad 0" component. A good fit is obtained with a peak at ~ 1600 MeV in the broad 
component, dropping smoothly to 2400 MeV, plus an interfering component over the mass range 
1850 to 2200 MeV, as shown below in Fig. 5(a). The smooth drop in the broad component from 
1600 to 2400 MeV is close to the P 3 dependence expected for point-like production. 

Although the magnitudes of the dominant components are well determined in Fig. 3, their 
relative phases are poorly determined. The final state with J p = 0~ is produced by an Ml 
transition and that with J p = + is produced by an El transition. The program takes the trace 
over the two photon polarisations allowed by Gauge Invariance and there is then no interference 
between 0~ and + . There is interference between + and 2 + and between 0~ and 2 + , but it is 
fairly small in practice. However, when the program takes the trace using all three spins 0~, + 
and 2 + , there does remain some significant information from details of angular correlations on 
the phase of the 0~ amplitude and this will be discussed in detail in Section 8. 

Fig. 4 displays information on phases from the slice fit after integrating over all angles. 
Errors for _ , + and 2 + are determined from the change in log likelihood as the phase of each 
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Figure 5: In each panel, the data points and upper histogram show the full An signal for 
reference. There is a cut from 1245 to 1325 MeV to remove /i(1285) — > An. Lower histograms 
show individual components: (a) Full 0~ intensity (full histogram) and from the broad CT alone 
(dotted); (b) 77(1970); (c) coherent sum of all + ; (d) coherent sum of all 2 + components (full 
histogram) and background (dotted). 



individual partial wave is moved in steps. Errors for the phase difference between 0~ and + 
span most of the range to 360°, so these results are inconclusive and are not shown. The 
relative phase between 2 + and + is shown in Fig. 4(a). The two 2 + amplitudes are compatible 
within errors with a constant or with a slowly rising phase difference in Fig. 4(b). However, 
Fig. 4 is not able to display the delicate correlations between all three spins if the angular 
information is used in addition. This point will be discussed in Section 8. 



5 Full fit to data 

The dominant feature in the 4n mass spectrum is the broad J PC = CT + signal in pp, extending 
over the whole mass range from 1500 to 2400 MeV. It will be described here as 'the broad 0~'. 
The rapid rise of the 0~ cross section at low mass is due to the threshold for pp. 

In Fig. 5(a), the upper histogram in each panel shows the overall fit to the 4-7T mass projection. 
The fit to the mass region 1500-1750 MeV is not perfect. There appears to be one low point at 
1610 MeV. No available ingredient is narrow enough to fit a single bin. The closest is / 2 (1565), 
but it has a full-width at half-maximum of 150 MeV. Section 8 will present improvements if the 
1600 MeV pp peak is fitted as a resonance rather than a threshold cusp. 

Fig. 5(a) shows as the lower full histogram the whole fitted 0" component. The dotted 
histogram shows the slowly varying 0" component, which falls smoothly from 1560 to 2400 
MeV. Fig. 5(b) shows the component fitted as a - resonance at 1970 ± 25(stat) ± 60(syst) 
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m(Jt + 7t"jc + JC ') (GeV) 

Figure 6: + contributions from (i) the coherent sum of /o(1370) + /o(1500), (ii) /o(1790) and 
(iii) / (1370). There is a cut from 1245 to 1325 MeV to remove /i(1285) -> 4vr. 



MeV with T = 210 ± 25(stat) ± 60(syst) MeV. Systematic errors cover the whole range of fits 
and variations in the form factors describing the broad 0~ component. The 77(1970) contribution 
is close to 90° out of phase with the broad 0~, producing a dispersive shaped curve in Fig. 5(a) 
due to destructive interference between the real part of the 77(1970) and the imaginary part of 
the broad 0~ amplitude. 

Fig. 5(c) shows the full + component, with clear peaks at 1500 and 1800 MeV. In Fig. 5(d), 
the full histogram shows the 2 + component, which may be fitted by a broad resonance at 2150 
MeV, plus a small contribution from / 2 (1270) — > pp. This panel also displays as the dashed 
histogram the total incoherent background, which is included in the fit to data. In more detail, 
the / 2 (2150) requires almost equal contributions from decays to pp and / 2 (1270)cr. The fitted 
mass is 2150 ± 29(stat) ± 60(syst) MeV and the width 506 ± 30(stat) ± WO(syst) MeV. The 
largest uncertainty in fitting parameters arises from the large and rapidly varying background 
of Fig. 5(d). 

The mass of the broad 2 + component is significantly higher than the / 2 (1940) found in BES 
I data [5] and also higher than the / 2 (1950) of the PDG, which quotes M = 1944 ± 12 MeV, 
T = 472 ± 18 MeV. This raises the possibility that the 2 + component observed here consists of 
two or more unresolved states. In the production process, the short-range interaction between 
the photon and gluons limits the process to production of 3 P 2 qq states and excludes 3 F 2 . 
Crystal Barrel analysis identifies two 3 P 2 states, namely / 2 (1910) of the Particle Data Group 
and / 2 (2240) [15], well identified by polarisation data for pp — > 7r + 7r~ , but incorrectly listed by 
the Particle Group under the ss state / 2 (2300). An alternative fit to the 2 + signal is made using 
/ 2 (1910) and / 2 (2240). This gives an almost identical fit with log likelihood better than the first 
alternative by only 0.8; this is not statistically significant. The fit to the An mass projection 
shows no significant difference either. So this is a feature of the data which cannot be resolved 
conclusively at the moment. The Table of results will use the single broad / 2 (2150). 

The fit to 0+ definitely requires contributions from all of / ( 1370), / (1500) and / (1790). The 
mass and width of / (1370) are fixed at values 1309 and 325 MeV from the recent re- analysis of 
Ref. [13]. If /o(1370) is omitted from the fit, log likelihood is worse by 100.3, a highly significant 
amount. Because of the rapid increase in 4n phase space, its peak in 4-7T is at 1390 MeV but it 
has a long tail in 4-7T at high mass and interferes with all other + components. Consequently, it 
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Figure 7: The poor fit to the An mass spectrum if the mass and width of / (1790) are changed 
to those of / (1710). 
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Figure 8: Fit to mass projection for (a) 2n, (b) all 3tt, (c) pn, with a cut on the p from 670 to 
870 MeV; (d) the azimuthal angle x between the planes of p pairs. 



is not well determined in magnitude. The possible process a — > An is inconsistent with the slice 
fit of Fig. 3(b) since it would give a broad contribution rising continuously with mass according 
to An phase space. 

Fig. 6 shows contributions from the coherent sum of /o(1370) and / (1500). The effect of 
/o(1370) is needed to broaden the contribution from / (1500). But the resulting tail interferes 
with /o(1790). If the parametrisation of the / (1370) component is varied within its known 
errors, the observed mass of the 1790 MeV peak and its full width at half-maximum remain very 
steady at M = 1800 ± 25 MeV, T = 230 ± 30 MeV. The observed peak is visually incompatible 
with /o(1710), to which the PDG assigns a width of 137 MeV. So the observed peak provides 
independent evidence for the existence of /o(1790). If the mass and width of /o(1790) are 
changed to PDG values for / (1710) MeV and the fit is re-optimised, the An mass projection is 
very badly fitted, see Fig. 7. 

An important result is that previous evidence for /o(2100) in these data has decreased sub- 
stantially; it is now almost invisible in the mass projection of Fig. 5(c). The reason for the 
change is the new evidence for the r/( 1970) . The / (2100) does appear clearly in E785 data on 
pp — > (7777)71"°, and there is independent evidence for it in Crystal Barrel data [15]. 

Fig. 8 shows further details of the fit to data. Fig. 8(a) shows the fit to the p(770) peak. 
The 37T mass projection of Fig. 8(b) and the pn mass projection of Fig. 8(c) exhibit the peak at 
1320 MeV due to the a2(1320)p background. Fig. 8(d) shows the azimuthal angular distribution 
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between and Kb-, requiring a large 0~ component, which varies as sin 2 \- Fig. 8(d) also 
reveals a substantial flat component. Unfortunately, the dependence on sin 2 x alone does not 
separate 0~, + and 2 + components cleanly. Note also that there are two p°p° combinations 
which interfere coherently; this further hampers a separation of the three available J p from sin 2 \ 
alone. This is a well known problem in analysis of multi-body final states. However, the full 
multi-dimensional correlations between production angles, decay angles of intermediate isobars 
and the azimuthal angle between pairs such as pp do separate J p cleanly, as displayed in Fig. 3. 
There is no simple way of displaying the 3-dimensional angular correlations (or 5-dimensional 
when two pp combinations are included); it is necessary to rely on the detailed amplitude analysis 
and the variation of log likelihood with parametrisations of individual resonances. 



6 More detail on formulae 

The final state X has a radius of interaction between final pp or 00 configuations of order 
R = 0.7 fm. It is necessary to include a form factor describing this radius of interaction. The 
conventional form factor for decays is 

F d = exp(-k 2 R 2 /6) = exp(-a d k 2 ), (1) 

where a d = (l/6)(i?/ftc) 2 = (l/6)(i?(/m)/0.197321) 2 with a d in (GeV/cT 2 . 

For decays with L = 1, the standard Blatt-Weisskopf centrigual barrier factor in the ampli- 
tude is 

B = K/[{hc/R') 2 + K 2 } 1 ' 2 . (2) 

Note that the derivation of this formula replaces the centrifugal barrier by an equivalent square 
barrier, so R' may be slighly different from R in the form factor. However, the Gaussian form 
factor is only an approximation to unknown wave functions, so it is convenient to fit data using 
a single radius parameter, setting R = R'. 

There is also a centrifugal barrier and form factor for the production process, except that R 
takes a much smaller value; the parameter a d of Eq. (1) is replaced by a much smaller parameter 
a p . The intensity of the J p = CT signal as a function of mass is given by the square of the form 
factor and centrifugual barrier factor, multiplied by a factor P for phase space in the production. 
The result is the factor P 3 quoted above, multiplied by a weak form factor F p = exp(—a p P 2 ). 

A general remark concerns the mass of the pp peak. If the width of the p were zero, the 
P-wave intensity for pp would rise from threshold as k 3 where k is the momentum of each particle 
in the pp system. A centrifugal barrier and decay form factor of conventional radius ~ 0.7 fm 
would make the pp and ujuj intensities both peak in the range 1700-1800 MeV. The width of the 
p spreads this peak out further than that for ujuj. It is the factor P 3 for the production cross 
section from a point-like source which reduces the mass of the peak to ~ 1600 MeV. 

It is necessary to fold the width of the two p into the evaluation of An phase space. This is 
given by 

y-h/i-W dsi f (VS-VH)* ds 2 8\k\ \h\ \k 2 \ im , 2 2 2 

PaA s ) = / \Ti{si)\ \T 2 {s 2 )\ exp(-2a d k ), (3) 

JAml 71 J4ml 71 y/S Si S 2 

where k\ and k 2 are momenta of pions from the decay of each p in its rest frame, and k stands 
for the momenta of the p in the rest frame of the resonance X. The amplitude T(s) is the usual 
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Figure 9: The dispersive term m(s) for the broad contribution. 



Breit-Wigner amplitude for the p. These integrals are done numerically and then parametrised 
by smooth functions of s. 

6.1 The effect of thresholds on resonances 

It is common practice to fit resonances with Breit-Wigner amplitudes with constant width. 
However, the Flatte formula shows this to be only an approximation. That formula describes 
the amplitude by 

T = F{s)/[M 2 -s-iY,9fpi(s)], (4) 

i 

where F(s) = F p F d , g { is the coupling constant and p; L the Lorentz invariant phase space for 
each channel %. 

Eqn. (4) however requires a further term if the amplitude is to be described correctly by an 
analytic function of s. It has been known since the 1950's that the full form for T is 

T = F( S )/[M 2 - S -m( S )-^^ 2 A ( S )] (5) 
mis) - ( g ~ M2 ) P f ^Pi{s')ds' 

m[s) ~ ^7 v-sW-Apy [ j 

The additional term m(s) describes the s-dependence of the real part of the amplitude going 
with the s-dependence of the imaginary part. It is essential in fitting the present data near the 
pp threshold; without it, the fit is very bad, see Fig. 11(b) below. The mass dependence of m(s) 
is shown in Fig. 9. The new parametrisation of the broad 0~ amplitude supercedes primitive 
alternatives used in Refs. [4] and [5]. 

A general comment is needed on the role of M 2 in Eqs. (4)-(6). There may be further decay 
channels, presently unknown. Also the exponential form factor of Eq. (1) is a conventional edu- 
cated guess; it is needed to make the dispersion integral of Eq. (6) converge. Realistically there 
are uncertainties in the high energy behaviour of this integral, well known as a renormalisation 
effect. The value of M 2 may be viewed as a way of applying an empirical correction to m(s) 
for such effects. A single broad resonance has been used in the earlier work by Bugg, Dong and 
Zou [20]. In the earlier work M was 2.19 GeV; now it moves to 2.04 GeV. 

Fig. 10 shows that similar threshold behaviour is observed in -qim , KK*, uouo (from Mark 
III) [16], K*K* and 00. Fig. 10(a) shows the broad 0~ component fitted to ^(rjTnr) data, the 
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sum of rja and aovr. Fig. 10(b) shows the coherent sum of 7/(1440) and the broad CT fitted 
to the KK* channel in an accompanying paper on J/\I/ — > 7(777r + 7r~) and / -f(K ± Kg7r T ) [17]. 
Fig. 10(c) shows what is observed with J p = 0~ in 7pp. A further threshold peak (not shown) 
is observed in BES II data for J p = + on J/$ ^ [18]. The natural explanation of 

all these processes is that J/^ — > 7GG, where G are gluons. These produce pairs of vector 
mesons via colour neutralisation. In turn, these may de-excite to lower mass configurations, e.g. 
[pp]L=i [w]l=i- This can explain why the r\o intensity peaks strongly at the same mass as 
pp. The inference is that all these peaks may be threshold cusps. Such cusps are well known 
in other processes. Examples are threshold peaks in nd — > NN and pp — > AA [19]. The other 
possibility is that one or more of these peaks is resonant. 

The data have been fitted using Eqs. (5) and (6) including the effects of cusps in all coupled 
channels. The coupling constants gf are adjusted iteratively to reproduce branching fractions 
of every channel in Fig. 10. As a result, the form of the broad 0" is tightly constrained. The 
dispersive term associated with the opening of the pp threshold plays a crucial role. Without it, 
the fit to data, shown on Fig. 11(b), is hopelessly bad. Fig. 11(a) shows the Argand diagram 
for Eq. (5). Note that what is plotted is the amplitude without the form factor F(s) of the 
numerator. The reason for this is that the factors P for production and k for decay are built 
into the tensor expressions for amplitudes, as well as the associated form factors. The intensity 
observed in data at low momentum in Fig. 10(b) is inflated by the factor P 3 for production. 
Note too that there are two interfering pp combinations in Fig. 11(b). 

Approaching 1900 MeV, the magnitude of the amplitude on Fig. 11(a) is decreased by 
the opening of the strong K*K* threshold. As this channel decreases above 2200 MeV, the 
amplitude recovers. Note that this threshold does not account for the structure near 1900 MeV 
which has been fitted by the additional 7/(1970). 

On the Argand diagram of Fig. 11(a), there is not a complete loop corresponding to a 
resonance in the mass range near 1600 MeV. It is being fitted as a simple threshold cusp. Below 
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Figure 11: (a) Argand diagram for the broad ; masses are shown in GeV; (b) the bad fit to 
the 0~ mass projection with the dispersive term m(s) omitted. 



the threshold, the real part of the amplitude moves positive. Over the threshold, the imaginary 
part of the amplitude increases rapidly. Above the threshold, the real part of the amplitude 
decreases slowly. The result is a half loop, characteristic of a cusp. 

However, from Eqs. (5) and (6) and fitted parameters, it turns out that there is a broad 
pole at M = M - iT/2 = 1.564 ± 0.012 - i(0.139 - i0.035) GeV; errors are systematic and are 
estimated from the parametrisation of the pp channel and its extrapolation into the complex 
s-plane. There is no discernable pole due to the weaker P-wave ujoj threshold which peaks at 
1750 MeV. There are more distant poles, including one at 1.862 — i0.210 GeV associated with 
the opening of the K*K* channel. 

For present data, the parameter controlling the form factor in Eq. (1) optimises at 
2.25 ± 0.25 GeV~ 2 , corresponding to a reasonable radius of interaction for decays, 0.73 ± 0.04 
fm. The value of a p for the production process optimises at 0.06 GeV -2 , corresponding to a 
radius of interaction 0.12 fm. As an aid to future work where the broad 0~ is relevant, the 
fully annotated Fortran code used for it is available directly from the author. This code includes 
numerical parametrisations of Pi(s) for ad in steps of 0.25 GeV" 2 . For values of ad or a p different 
to those used here, it will be necessary to re-optimise gf fitted to all relevant sets of data. 

7 Branching fractions 

Table 1 summarises the significance levels for each component, measured by changes AS in 
log likelihood when the component is removed from the fit and all other components are re- 
optimised. Statistically, x 2 is approximarely twice AS. However, experience shows that signals 
less than 20 in log likelihood are of questionable significance because of the possibility of sys- 
tematic errors. Contributions in this category are the weak decays of / (1790) and / (2100) to 
pp and also / 2 (1565), which produces a change of log likelihood of just 11.9. 

The total number of selected events is 23591. The total background fitted up to a mass of 
2.4 GeV is 27.3% of all events, leaving a total signal of 72.7%. The corresponding branching 
fraction for all events is 

J/V -> 7(7T + 7r"7T + 7r") = (25.5 ± 2.6(syst)) x 10~ 4 , M(4tt) < 2.4 GeV. (7) 

The error is systematic and arises mostly from uncertainty in the experimental background. In 
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channel 



5{lnL) Branching fraction xlO' 



77(2190) -> pp 
?7(1970) -> pp 



> 1000 

233.5 2.0 ±0.3 

13.7 ± 1.4 



Both 0" 



/o(1370) -> aa 
/o(1500) -> ad 
/o(1790) -> ffff 
/o(1790) - pp 
/o(2105) -> ffff 
/o(2105) - pp 
/ 2 (1565) - pp 
/ 2 (2150) - pp 



100.3 1.1 ±0.4 

163.8 1.5 ±0.4 

140.0 5.8 ±1.2 
3.0 



25.2 
6.0 



11.9 0.7 ±0.3 

70.2 



/ 2 (2150) -> / 2 (1270)<7 
/ 2 (2150) ^pp ±/ 2 (1270) 



(7 



178.5 

190.5 3.1 ±0.5 



Table 1: Changes in log likelihood when each channel is dropped from the fit and remaining 
contributions are re-optimised. Column 3 shows branching fractions for J/\f r — > 7X, X — > 
7r + 7r _ 7r + 7r _ . Missing elements are not quoted because of large uncertainties due to interferences. 
Errors include systematic errors on the overall normalisation. 



determining branching fractions given in Table 1, interferences within one channel between the 
two possible 7r + 7r~ combinations are included, but interferences between different channels are 
dropped. This is because interferences within one resonance should be universal, but interfer- 
ences with other channels vary from process to process. Note that the PDG quotes branching 
fractions for all pp and 00 charge states. These are conventionally obtained by multiplying the 
branching fractions of Table 1 by a factor 3 for pp and a factor 9/4 for 00. However, it is likely 
that interferences between different charge configurations of a single resonance will be different 
between p°p° and p + p~ . Chen et al. [21] have estimated the magnitudes of such effects. 

The fitting program itemises the magnitudes of all individual resonances and all interferences 
between them. A warning is that interferences are quite large within one J p '. There is destructive 
interference between the two decay modes of / 2 (2150), so the branching fraction is quoted only for 
their coherent sum. There is also destructive interference between the broad 0~ amplitude and 
77(1970), so the table quotes the total 0~ branching ratio. Within the + sector, the branching 
fractions of /o(1500) and /o(1770) are quite well determined, but there is a poorly determined 
interference with the tail of /o(1370). The magnitude of any possible /o(2100) contribution is 
difficult to determine because of interferences with the high mass tails of all of / (1370), / (1500) 
and / (1790), and is therefore not quoted. 

There is no direct interference between J p = 0~ and + . Interferences between + and 2 + 
are also quite small. The summed events for 0~, + and 2 + add up to 17166 events, to be 
compared with the total of 17146 events after background subtraction. 



15 



I 



1200 
1000 
800 
600 
400 
200 




1200 
1000 
800 
600 
400 
200 










L (a) 


: ¥ ++ 










1 


.5 2 



m(7l + 7t"7t + 7t") (GeV) 









L (C) 


: +F 








^ ■ 


1.5 2 



m(n + KK + n) (GeV) 



1200 
1000 
800 
600 
400 
200 




1200 
1000 
800 
600 
400 
200 










L (b) 


r+4i 

++ 










1 1 


.5 2 



m( 7t 7C") (GeV) 









L (d) 














k^4^ ^ 




1 


.5 2 



m(7t + 7l"7l + 7t") (GeV) 



Figure 12: The fit, including an 77(1560), to the An mass projection and magnitudes of fitted 
components in the same format as Fig. 5. 



8 Subtle evidence for a + resonance near 1560 MeV 

The 4ir mass projection shown on Fig. 5 is not fitted perfectly from 1500 to 1650 MeV. There is 
one low point at 1610 MeV. Moving the edges of bins or changing the bin width does not alter 
this point significantly. There is also some evidence for a similar dip in DM2 data [2] and a dip 
in Mark III data at 1650 MeV [4]. 

It is unreasonable to invent a new narrow resonance to fit the single low point at 1610 MeV. 
However, the observed pole in the broad 0" component is already a hint that the pp threshold 
peak is resonant. This is what would be called in today's terminology a dynamically driven 
pp resonance. It is quite likely to mix with the nn radial excitation expected in this general 
mass range, for example the radial excitation of 77(1295) if it exists. Suppose there is such a 0~ 
resonance in the range 1550-1700 MeV. Experience with other sets of data, for example Crystal 
Barrel data, Refs. [13] and [22], is that a resonance at or close to a broad threshold is in practice 
described well by a simple pole with a high mass dispersive tail given by the term m(s) of Eq. 
(6). That is, a good approximation can be obtained by adding a simple Breit-Wigner pole to 
the dispersive amplitude used so far. 

It is instructive to fit the data this way, scanning the mass and width of the additional pole. 
It turns out that this gives a distinct improvement in the fit to the 4-ir mass projection and the 
dip at 1610 MeV. It gives an improvement in 2x log likelihood (which follows a distribution 
very close to x 2 ) °f 194 for four extra fitting parameters, statistically a 10.6 standard deviation 
effect. Systematic errors quoted for mass and width cover all variations as masses and widths of 
/o(1500) and /o(1790) are varied over their allowed ranges from earlier data; systematic errors 
also cover uncertainty in the dispersive contribution of the rjirn channel. 
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The improvement in the fit comes partially from the fit to the An mass spectrum. However, 
it also comes from systematic improvements amongst decays to + , 2 + and CT partial waves. 
The fitting program takes the trace over both allowed spin states of the photon. There are 
interferences between 0~ and 2 + and between + and 2 + , dependent in subtle ways on three- 
dimensional angular correlations, hence giving information on relative phases. 

The new fit is shown in Fig. 12 in precisely the same format as Fig. 5. There is an obvious 
improvement in the fit to the 4-7T mass spectrum from 1500 to 1650 MeV, particularly near the 
dip. The signal due to / (1500) in Fig. 12(c) decreases slightly and the dip before / (1790) is also 
slightly lower. No change is required to the parameters of / (1790). The mass of 77(1970) does 
not change, but its width decreases from 210 ± 25(stat) ± 60(syst) to 165 ± 25(stat) ± 60(syst) 
MeV. A knock-on effect of this reduction is that the width of the broad 2 + component at 2150 
MeV decreases from 506 MeV to 300 ± 75 MeV. The evidence for / (2100) increases slightly 
on Fig. 12(c). A definite / 2 (1565) appears in Fig. 12(d). Except for this state, branching 
fractions remain within the errors shown in Table 1. The branching fraction for /2(1565) rises 
to (1.3 ± 0.3) x KT 4 , nearly double that of Table I. 
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Figure 13: The Argand diagram for the coherent sum of the 1560 MeV 0~ pole and the dispersive 
cusp; masses are shown in GeV 

The Argand diagram for the coherent sum of the 1560 MeV pole and the dispersive cusp is 
instructive and is shown in Fig. 13. The pole term is ~ 35° out of phase with the dispersive 
term, and the result is to increase the speed with which the amplitude describes the Argand 
loop of Fig. 13. As one sees from Fig. 12(a), the direct effect on the whole 0~ amplitude is to 
produce half the required dip at 1610 MeV. The remaining dip is created by a decrease in the 
+ amplitude between / (1500) and / (1790). The increased / 2 (1565) intensity in Fig. 12(d) 
and its interference with the + amplitude contributes to the fit to the 4n mass spectrum near 
1560 MeV. 

The coherent sum of the additional pole and the broad 0~ has been parametrised as a function 
of s and extrapolated into the complex plane. A pole is observed at M — 1560 ± 12 — i(143±35) 
MeV, close to that of the broad 0" alone. 

The improvement in log likelihood when the 0~ pole at 1560 MeV is included is larger than 
for most of the other contributions in Table 1, so the evidence for this resonance must be taken 
seriously. In summary, there is no doubt about the existence of a strong pp peak at ~ 1600 
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MeV. However, evidence for resonant behaviour is hidden in the angular correlations. It would 
be valuable to have further 4ir data in which the 0~ component interferes with another well 
defined component; central production in the Compass experiment is one such possibility. 

9 A conflict with jujuj analysis 

A BES II paper on J/$ ^ ^uouo claims a large and almost pure 7/(1760) signal in uouo [23]. No 
such resonance is observed in the present data. This raises a problem. There is a well known 
relation, coming from SU(2) symmetry, that an I = resonance should have equal couplings to 
uouo and p°p°. Physically, light quarks do not discriminate between charges and therefore couple 
equally strongly to uouo and p°p°. Because there are three charge states for pp, the relation 
is normally written g 2 (pp) = 3g 2 (uouo). This relation applies equally to qq states, hybrids and 
glueballs, since it springs from SU(2) symmetry, which is obeyed by all these configurations. The 
Appendix gives an elementary derivation of this result for a qq state. Similar algebra generalises 
the result to all three types of state. 

The branching fraction quoted for production of 7/(1760) in the ^uouo data is (1.98±0.08(siai)± 
0. 32 (syst)) x 10~ 3 . This is larger than the branching fraction in present data for production of 
p°p° with J p = 0- over the entire mass range: (1.37 ± 0.14) x 10 3 . It should lead to a huge 
peak with resonant phase variation interfering with the broad 0~ component. The present data 
are in complete disagreement with that possibility. Furthermore, such a large resonance is also 
in conflict with all earlier data from BES 1 [5], Mark III [1] and DM2 [2]. The DM2 group 
did claim a small 7/(1760) signal originating from the peak near this mass in the total An mass 
projection; however, this had a branching fraction an order of magnitude smaller than BES II 
claim. Furthermore, DM2 did not try a + — > oo contribution which fits the peak in present 
work. 

The likely explanation of the conflicts is that the 0~ signal in uouo is not a resonance, but 
a threshold cusp. Note that there is no phase information identifying resonant behaviour of 
7/(1760) in BES data. It was simply assumed that the peak is resonant. A threshold cusp arises 
because the uouo intensity rises initially from threshold as k 3 , but peaks near 1700 MeV due to 
the form factor exp — (o^p 2 ); here k is the momentum of each uo in the uouo rest frame. The 
precise position of the peak depends rather strongly on eta- The BES publication does not give 
details of the form factor assumed for the P-wave. 

This explanation may be tested by refitting the uouo mass projection of BES II data using the 
uouo decay of the broad 0~ fitted to present data. The 7/(1970) component should also obey the 
SU(2) relation between g 2 ^ and g 2 o p o- This constraint is applied by fixing relative magnitudes 
of 7/(1970) and the broad 0~ to the value fitted to present data. However, the phase of 7/(1970) 
relative to the broad - is fitted freely, since it is determined by multiple scattering amongst all 
components, which are different in uouo and 4-7T final states. 

Decays of /o(1790) to uouo are also included, but are limited in intensity to 20% of the broad 
0~. This is the maximum which can be fitted to / (1790) — > p°p° in present data. Decays of 
the + state produce a cos 2 x decay distribution (where x is the angle between uo decay planes) 
which combines with an equal amount of 0~ to produce a flat component in the x distribution. 
This allows the 0~ component of BES II data to rise to ~ 46% of those data. 

There may be small 2 + contributions from / 2 (1565) and / 2 (1910)//2(2150). Figs. 5(d) and 
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Figure 14: The fit to BES II ^ujuj data as a threshold cusp. In (b), the upper full curve shows 
the total contribution; the dotted curve shows the broad 0~ alone, and the dashed curve its 
coherent sum with 77(1970). The chain curve shows the / (1790) intensity and the lower full 
curve the experimental background scaled from the BES II publication. 



12(d) show that these are almost zero at 1750 MeV. They can only be determined from the ^ujoo 
data using the full angular correlations. They are therefore omitted from the test made here. 

Fig. 14(a) shows the resulting fit to the uu mass spectrum. The optimum fit requires some 
reduction of the parameter of the form factor for decays to 2.0 GeV -2 . The fit is close to the 
data, showing that they may indeed be fitted by a non- resonant threshold cusp. The optimum 
fit would require ad < 2.0 GeV -2 , but that would create an inconsistency with the pp data. 
Details of the components are shown on Fig. 14(b). The full curve is the same as in Fig. 14(a) 
as a reference. The dotted curve shows the contribution from the broad 0~, and the dashed 
curve shows its coherent sum with the 77(1970). The overall effect of r/(1970) is quite small. The 
chain curve shows the /o(1790) contribution, fitted with the required phase space for S-wave 
decays. The broad - makes a small but significant contribution above 2 GeV; BES II fitted 
that mass range purely as background. To allow for this, the s-dependence of the background 
is taken from the BES publication, but its absolute magnitude is fitted freely and is shown by 
the lower full curve on Fig. 14(b). 

The fit in Fig. 14 is made without considering the absolute magnitude of the branching 
fraction for J/ty — + ^ujuj. The BES II publication quotes branching fractions in their Table 2 
adding to (2.77±0.12) x 10 -3 of all J/^f decays. It is a matter of conern that this is considerably 
larger than quoted by DM2 [2], 1.42 ± 0.2 ± 0.42 and Mark III, 1.76 ± 0.09 ± 0.45 [1]. Also, 
one would expect the branching fraction for pp after integration over all masses to be roughly 
the same as for uu; the pp peak is simply spread out compared with uuo by the width of the 
p. The disagreements over branching fractions raises the question of backgrounds in the BES 
data. Those used in the BES publication were estimated by Monte Carlo simulations, which 
really need checking against the standard sideband subtraction procedure. Also the BES II 
analysis chose the best of 15 combinations of photons fitting ^ujuj; wrong combinations will 
make a significant background, but this background is not discussed in the BES publication. 
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10 Summary and Conclusions 



In summary, the new data confirm earlier detection of + states at 1500 and 1790 MeV in this 
channel. There is a broad 2 + signal at ~ 2150 MeV decaying to both pp and / 2 (1270)<7, but 
this may be an unresolved combination of more than one state. 

For J p = 0~ the new data give a distinctly improved definition of the broad component. 
The exponential form factor for decays is well determined, with aid = 2.25 ± 0.25 (GeV/c) 2 , 
corresponding to a radius of interaction of 0.73 ± 0.04 fm for decays. The exponential form 
factor for production is also well determined, with very small a p = 0.06, corresponding to a 
radius of interaction of 0.12 fm for the production process; this is consistent with the expected 
point-like interaction. There is also a definite dispersive shaped effect with a resonant phase 
variation with respect to the broad - component, requiring the presence of a 0~ resonance 
with mass 1970 ± 25(stat) ± 60(syst) MeV and width 210 ± 25(stat) ± QO(syst) MeV, in fair 
agreement with a less definitive Crystal Barrel observation [15]. 

There is subtle but persuasive evidence that the observed pp threshold peak with J p = 0~ 
resonates at ~ 1560 MeV. Because there is only limited phase information from interference 
with / 2 (1565) (which in turn interferes with / (1500) and /o(1790)), confirmation of a resonance 
is needed in other data where some well defined component plays the role of an interferometer. 

In both Crystal Barrel data on (a) pp — > r]iT iT [24] and (b) pp — > 3n [25] , there is evidence 
for a J PC = _+ resonance in / (1500)t7 at 2285 ± 20 MeV in (a) and at 2320 ± 15 MeV in 
(b). The latter is a particularly simple channel exhibiting a clear peak. One can now identify a 
Regge trajectory of conventional slope for 7/(550), 7/(1295), 7/(1970) and 7/(2320). The possible 
7/(1560) falls 100 MeV below this trajectory. However, its mass is subject to a systematic error 
of ±60 MeV at present. 

I wish to acknowledge financial support from the Royal Society and Queen Mary College. I 
am also grateful for help from members of the BES collaboration in processing data and running 
the Monte Carlo simulation of acceptance and backgrounds. 

11 Appendix 

The 1 = nn combination is (uu + dd) / \[2. Its decay is via the production of a 3 Pq pair 
(uu + dd)/\/2. The final states is (uu + dd)(uu + dd)/2. An tutu pair has this compositon by 
inspection. 

What about pp decays? The isospin Clebsch-Gordan coefficients for an / = resonance 
coupling to two 1=1 particles is (p + p~ + p~ p + — p°p°) / a/3, where p + = ud, p° = (uu — dd) / \[2 
and p~ = —ud. So 

(p + p~ — p° p° + p~ p + )/V3 = — udud — (uu — dd)(uu — dd) /2 — udud (8) 

= (—2udud — uuuu — dddd + udud + udud — 2udud)/2 (9) 
= — (uuuu + dddd + udud + udud)/ '2 (10) 
= -(uu + dd)(du + dd)/2. (11) 

So the decay of / = nn produces the combination (p~p° — p + p~ — p~p + ), i.e. the same pp° 
intensity as ujuj and also p~ p + ■ 
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